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Abstract 

We propose that the stability of dark matter is ensured by a discrete subgroup 
of the U(1)b-l gauge symmetry, Z2(B — L). We introduce a set of chiral fermions 
charged under the U(1)b-l in addition to the right-handed neutrinos, and require 
the anomaly-cancellation conditions associated with the U(1)b-l gauge symmetry. 
We find that the possible number of fermions and their charges are tightly con- 
strained, and that non-trivial solutions appear when at least five additional chiral 
fermions are introduced. The Fermat theorem in the number theory plays an im- 
portant role in this argument. Focusing on one of the solutions, we show that 
there is indeed a good candidate for dark matter, whose stability is guaranteed by 
^2(B-L). 



The presence of dark matter (DM) in the Universe is clear observational evidence for 
the need of physics beyond the standard model (SmJ]. The DM abundance has been 
measured with a very high accuracy as [2] 

VtBMh^ = 0.1109 ±0.0056, (1) 

where h is the present Hubble parameter in units of lOOkm/s/Mpc. If the DM is made 
of as-yet-undiscovered particles, they must be electrically neutral and long-lived. In par- 
ticular, their lifetime must be much longer than the present age of the Universe, and the 
longevity may be ensured by a symmetry such as a Z2 symmetrjj^. There are, however, no 
exact global continuous or discrete symmetries, according to the argument on the quan- 
tum gravity by Banks and Seiberg [4J. Therefore, it may well be that the symmetry is an 
unbroken subgroup of a gauge symmetry. 

The U(1)b-l gauge symmetry is a very attractive symmetry beyond the SM, since 
it predicts three families of right-handed neutrinos. The observed small masses of light 
neutrinos are naturally explained by the seesaw mechanism [5], provided that the B-L is 
broken at a very high energy such as the grand unified theory (GUT) scale. Interestingly, 
a Z2 subgroup of the U(1)b-l gauge symmetry remains unbroken in the low energy, if the 
U(1)b-l is spontaneously broken by a Higgs field of charge 2, coupled to the right-handed 
neutrinos to generate the large Majorana masses. Here we define the normalization of 
the B-L charge so that the right-handed neutrinos carry the charge —1. In this letter we 
consider a possibility that the ^2(6 — L) is responsible for the stability of DM. 

The anomaly-free conditions associated with the U(1)b-l symmetry are satisfied for 
the SM fermions plus the three right-handed neutrinos, which however contain no DM 
candidatqj- Accordingly, we need to introduce an extra particle, which must be electrically 
neutral and (quasi-) stable to be a good candidate for DM. In this letter we consider a 
set of SM-gauge- singlet chiral fermions charged under the U(1)b-l, some of which are 
to be the DM. One important constraint on those fermions is that they must satisfy the 



^ One exception is a primordial black hole which is the only candidate for DM in the SM. 
^ Another possibility is a composite dark matter scenario [3]. 

^ If one of the right-handed neutrinos is extremely light, it can account for DM. In fact, it is possible 
to split the mass scales of the right-handed neutrinos while keeping the success of the seesaw mechanism 
in the the extra dimensional framework [B] . 
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anomaly-free conditions associated with the U(1)b-l, which severely limit the possible 
number of fermions as well as their (relative) B-L charges. We assume the B-L charges 
to be rational, since otherwise we would have an exact global symmetry in contradiction 
with the quantum gravity In particular, by observing that all the SM fermions and 
the right-handed neutrinos have odd B-L charges, the lightest new fermion with even 
B-L cha.,e wU. be .tab.e and the^efo. a cand.date fo. DmS Also, if the B-L charge is 
relatively large, the fermion mass tends to be light enough to be produced in the early 
Universe. 

In this letter, we show that one of the "minimal sets" of new chiral fermions contains 
those with large B-L charges. Interestingly, one of them turns out to be a good candidate 
for the DM in the Universe. Here we define the minimal set of fermions such that it 
contains no two different chiral fermions which have the same B-L charge: i.e. the number 
of the Qi charged chiral fermions ipi is one or zero. Thus, this is one plausible example in 
which the stability of DM is ensured by a discrete subgroup of a local symmetryl^ 

Let us first see how the anomaly-free conditions limit the possible number and charges 
of the additional fermionslll There are two anomaly-free conditions: one is from [[/(1)b~l]'^ 
anomaly and the other from the gravitational [f/(l)B-L] x [graviton\^ anomaly. They are 
given by 

EQHO, (2) 

i=l 
n 

i=l 

where we assume Qi 0, and n is the number of newly introduced fermions, ipi [i = 1- 
n). We define all the new chiral fermions to be right-handed fermions without loss of 
generality. According to the argument of Banks and Seiberg [3], Qi must be a rational 



''The B-L corresponds to the fiveness U(l) in the SU(5) GUT. Notice that all the SU(5)-invariant 
operators of quarks and leptons (5*, 10) carry integer charges of the fiveness U(l). 

^ This applies to the case that the B-L charges of the additional new fermions are integer. 

^ Other models can be found in Ref. [7] , in which the stability of DM is ensured by a discrete symmetry, 
a remnant of hidden U(l) gauge symmetry. 

Mathematical analysis for finding a set of chiral fermions satisfying the anomaly-free conditions was 
given in Ref. |8j. 
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number. By multiplying a non-zero integer a with the both sides of Eqs. ([2]) and ([3]), 
therefore, we can always rewrite them in the following form, 

n 

Y.^Z,f = 0' (4) 

i=l 

n 

i=l 

where {Zi, ■ ■ ■ Z^} = {aQi, ■ ■ ■ , aQ„} are integers. For convenience, we rearrange {Zi} 
in ascending order: Zi < Z2 < ■ ■ ■ < Zn. Let us seek solutions of Eqs. (jl]) and with 
n > 1. The solutions of Eqs. ([2]) and ([3]) can be obtained by multiplying the integer 
solutions by some rational number. 

In the case of two new fermions {n = 2), one can easily find that the solutions of 
Eqs. (jl]) and ([5]) are given by Zi = — Z2. In this case we can construct a mass term, 

C = ^mipiip2 + h.c. (6) 

The mass parameter m is not constrained by any symmetry and hence can be as large as 
the Planck scale. This does not lead to interesting phenomenology within our framework, 
and hence we do not consider the case of n = 2. In the following we will exclude vector-like 
fermions from the solutions. 

In the case of three new fermions {n = 3), we can easily see that there is no solution 
to Eq. (jll), by noting the famous Fermat theorem in the number theory [9j@ This fact 
requires at least four new fermions to cancel the anomalies. Thus we already suspect that 
such a solution, if it exists, may contain chiral fermions of large B-L charges. 

In the case of four new fermions (ra = 4), it is easy to prove there is no phenomeno- 
logically interesting solution to the both conditions, (j4]) and ([5]). If we erase Z4 from the 
two equations, we find 

{Z^ + Z2){Z2 + Z,){Z■, + Z,) = 0. (7) 

* Euler discovered a proof of the Fermat theorem in the special case of the cubic equation (i.e. Eq. ([4]) 
with n = 3) |10j . which however contained an important missing step. The complete proof was given by 
Kausler [TT] and many others. 
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-25 


-8 
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18 
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Table 1: Independent solutions to Eqs. (jlj) and ([5]) for max{|Zj|} < 25 for n = 5. 



Thus, two of t/'i, '\\)2 and ■i/'s must have the charges of opposite sign and equal magnitude, 
and so do the other fermion and ^/'4. The solution is therefore given by two pairs of 
vector-like fermions, in which we are not interested, as in the case of n = 2. 

Thus we are led to consider the case of at least five new chiral fermions in > 5). 
Introducing five additional fermions {n = 5) is the minimal non-trivial extension. The 
independent solutions to Eqs. (jll) and ([5]) with n = 5 and max{|Zi|} < 25 are given in 
Table [H We find that even this minimal extension leads to relatively large B-L charges. 
We show in the following that the minimal charge solution in Table [1] contains a good 
candidate for the DM. 

The seesaw mechanism [5] for neutrino mass generation suggests the Majorana mass 
of the (heaviest) right-handed neutrino at about the GUT scale. For this purpose, we 
introduce a Higgs field $ with the B-L charge 2. We assume that $ develops a vacuum 
expectation value (vev), ($) = vb-l, where vb-l represents the B-L breaking scale. The 
three right-handed neutrinos, A^^i, N2 and A^3, acquire a mass from 

^K3<^N3N3 + ^K2<!>N2N2 + ^Ki<l>NiNi, (8) 

where ^1,2,3 are coupling constants with \k,s\ > I ^2! > l^^il, and we have adopted a basis 
such that the right-handed neutrino mass matrix is diagonalized. The breaking scale can 
be inferred from the neutrino oscillation data, assuming that the couplings of the heaviest 
right-handed neutrino are of order unity. Then the B-L breaking scale vb-l is estimated 
to be about lO^^GeV. As noted before, there is unbroken Z2(B — L), as long as the 
U(1)b-l symmetry is broken only by the vev of $. Since Z2(B — L) is a subgroup of the 
gauge symmetry, it can be an exact symmetry which ensures the stability of DM. 

First, let us focus on the integer solutions to Eqs. ([2]) and ([3]), which are obtained by 
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identifying with ±Qi, i.e. by setting a = ±1. To be explicit we consider one of the 
integer solutions, 



(Qi, Q2, Q3, Q4, Q5) = (-9, -5, -1, 7, 8), (9) 

which corresponds to the first solution in Tabled] with a = 1. In this example, the chiral 
fermion with a charge 8, ip^, is stable because this is the only chiral fermion with an even 
B-L charge and it has no mixings with the other fermions. In other words, the stability 
of ip5 is guaranteed by the ^2(6 — L). Therefore the fermion ip^ is a prime candidate for 
DM in this solution. 

The ip5 acquires a mass from the following non-renormalizable operator, 

^^5^5, (10) 

where M is a cut-off scale. Then the mass of ip^ is given by 

where we have set the cut-off scale M to be the Planck mass, Mp ^ 2.4 x 10^^ GeV. In 
the following we will take f b-l = 3 x 10^^ GeV and M = Mp as reference values unless 
otherwise stated. The mass of ip^ is sensitive to the B-L breaking scale, but it can be 
heavy enough to be consistent with the Lyman-alpha data [12]. 

In order to account for the DM, ipr, must have the correct abundance. The main 
production process is pair production of ip5 from the SM fermions in plasma through 
the s-channel exchange of the B-L gauge boson [6j. The production is most efficient at 
reheating. The iJj^ number to entropy ratio can be roughly estimated as 



45 



T=Tr 



where H is the Hubble parameter, counts the relativistic degrees of freedom at the 
reheating, {av) ~ T"^ h^-h the production cross section, rij ~ is the number density 
of the SM fermions in plasma, Tr is the reheating temperature and the first equality is 
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evaluated at the reheating. The numerical solution of the Boltzmann equation gives a 
consistent result [13j. The DM abundance is given by 

Thus, the reheating temperature as high as 0(10^^) GeV is needed to account for the DM 
density by this production process. The thermal leptogenesis [H] works with such a high 
temperature. 

Let us comment on the properties of the other fermions, ipi^ip2-,'4''ii ^ind '?/'4, which have 
charges of —9, —5, —1 and 7, respectively. We fist note that will have a Majorana 

mass of order v-q-l similarly to the right-handed neutrinos. We can also see that ip^ and 
one linear combination of ipi and ip2 acquire a mass of order ^^b-l from the mass terms, 
$V'i''/'4 $*-?/'2^4; these heavy modes do not significantly contribute to the light neutrino 
masses, and so they can be safely integrated out. Also they do not have any cosmological 
effects, since they are not produced in the early Universe if the reheating temperature is 
smaller than their masses. On the other hand, the other orthogonal combination of ipi 
and ip2i denoted by ipi in the following, remains light. The ipi obtains a mass of order 
v^_i^/M^ ~ 7TeV both from the Majorana mass term and from the mixing with the 
right-handed neutrinos. The ipn quickly decays into a lepton and a Higgs boson before 
the big bang nucleosynthesis, and therefore has no drastic effects on cosmology]^ 

Let us comment on a possible problem with other solutions. For instance, we consider 
the first solution in Table[T]with a = —1. Similarly, we would then have a light fermion ip'^. 
Problem is that ip'^ would generally have larger mixings with the right-handed neutrinos, 
which spoils the seesaw formula. Thus, for such solutions, dangerous mixings between 
ip'f, and the right-handed neutrinos must be suppressed by introducing an approximate 
matter parity or extra dimensional set-up, in order to retain the seesaw mechanism. 

Next let us consider a rational number solution to Eqs. (|2]) and (|3]). One advantage 
of rational number solutions is that possible dangerous mixings with the right-handed 
neutrinos can be suppressed. On the other hand, the stability of the lightest Z2(B — L) 
even fermion is not necessarily ensured because of the presence of fermions of fractional 

^ The lepton asymmetry produced by the right-handed neutrino decay is not washed out if the mixing 
of tpi with the right-handed neutrinos is suppressed by 0(0. 1). 
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charge. Let us consider one example. Multiplying the previous solution with 1/2, we 
obtain 

(Qi, Q2, Qs, Q4, Qs) = -| ^, 4^ • (14) 

As before, ip^ is a candidate for DM and it has a mass of v^_j^/M^ ~ 6 x 10^ GeV for the 
same reference values of f b-l and M as those in the previous case. The correct abundance 
is realized if the reheating temperature is Tr ~ 8 x 10^ GeVo The thermal leptogenesis 
works for such high Tr. The other additional fermions have fractional charges and we 
can see that they remain massless and their contributions to the extra radiation will be 
negligibly small. Moreover, since they have fractional charges, they do not mix with 
the SM fermion and the right-handed neutrinos. We note here that ip;^ is not absolutely 
stable, because of the presence of fermions of fractional charge. In fact it mainly decays 
into ip2, V^a, the SM lepton and the Higgs boson with a lifetime more than twenty orders 
of magnitude longer than the present age of the Universe, for the reference values of f b-l 
and M. The hfetime is proportional to high powers of vb-l, and so, the decay could 
contribute to the high energy cosmic-ray spectrum for a slightly large value of wb-l- Note 
that another solution can be obtained by multiplying ( IT4l) with an odd integer, which 
however leads to too light DM mass. 

Finally we briefly discuss other solutions. Let us focus on the integer solution of 
{Qi,Q2,Q3,Qa,Q5) = (—18,-17,1,14,20). There are three ^2(6 — L) even fermions, 
one of which may be a candidate for DM. For a similar reason discussed above, one linear 
combination of ip4^ and ip5 remains light. Its mass is however of the order of ~ v^_-^/M^^ , 
and it is extremely light. The ^'2(6 — L) odd fermion, ip2, is also extremely light and 
it can be regarded as a stable particle in a cosmological time scale. Thus the light and 
stable particles in this solution are too light to be DM. Similar consideration leads us to 
conclude that additional fermions tend to be extremely light and DM candidates do not 
likely exist as the U(1)b-l charges become larger. On the other hand, these light stable 
particles may constitute a part of extra radiation in the Universe. If the number of such 
light fermions is large enough, they may be able to explain the extra radiation of the 



Since the number-to-entropy ratio is proportional to WbIl (^^^ (fT2l)). the DM rehc density ([T3|) does 
not depend on "yB-L- Thus the ips mass can be lighter if we choose a smaller value of vb-l- 
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Universe indicated by recent observations [21 [161 US HE] ■ 

Let us comment on other implications of our scenario. So far we have imphcitly 
assumed that the U(1)b-l symmetry is spontaneously broken during inflation and it is 
not restored after inflation. However, it is possible that the U(1)b-l gets spontaneously 
broken after inflation, if the mass of $ is light enough. Then the cosmic strings will be 
formed at the phase transition, which can be probed by the CMB measurement such 
as the Planck satellite. In general the DM abundance is modified in this case, but the 
estimate (fT2|) remains valid as long as the phase transition occurs at T > 10^^ GeV. 
Another implication of the integer solution (|9]) is that the DM mass is in the keV range 
so that it can affect the structure formation at small scales. We note that the lifetime 
of DM depends on the assumption that U(1)b-l is spontaneously broken only by $. If 
^2(3 — L) is also broken by another Higgs field of a B-L odd charge, the DM decays into 
the SM particle, which may produce observable signature in the indirect DM search. For 
instance, the DM decay produces an X-ray photon in the case of ([9]), while high-energy 
gamma-ray, anti-protons, and positrons are produced in the case of f|T^ . In fact, even 
without the additional Higgs, the DM is not absolutely stable in case of dHD, and the 
decay can contribute the cosmic rays if the B-L breaking scale is several times larger 
than our reference value. On the other hand, no signals are expected in the direct DM 
search. We have seen that a high reheating temperature is needed in order to obtain a 
correct DM abundance (I12p in the case of Since the inflaton mass must be larger than 
the reheating temperature [20], the chaotic inflation [2T| is one of the possibilities. This 
may be confirmed by the observation of the B-mode in the cosmic microwave background 
polarization. 

In this letter we have considered a possibility that the longevity of DM is guaranteed 
by the Z2(B — L). To this end we have introduced n SM-gauge-singlet chiral fermions 
charged under the U(1)b-L5 in addition to the three right-handed neutrinos. We have 
shown that more than four additional chiral fermions are required for the cancellation of 
the anomalies associated with the U(1)b-l- In the case of = 2, two fermions must have 
charges of opposite sign and same magnitude, and they can be integrated out if their 
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mass is as heavy as the Planck scale. Importantly, the Fermat theorem excludes the case 
of = 3, which have forced us to consider n > 3. In the case of = 4, the four fermions 
are divided into two pairs of vector-like fermions, which have charges of opposite sign and 
same magnitude, and this case was not of our interest. Thus, we are led to consider more 
than four additional chiral fermions. Interestingly, one of the minimal set of the integer 
B-L charges satisfying the anomaly-free conditions contains a good candidate for DM. The 
DM abundance can be explained if the reheating temperature is as high as 0(10^^) GeV 
or O(IO^) GeV, depending on whether the charges are integer or fractional. Both cases 
are consistent with the thermal leptogenesis scenario. 
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